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, Abstract. We consider the scaling behavior of directed percolation and of 

£f*) • the pair contact process with a conjugated field. In particular we determine 

f— *^ ' numerically the equation of state and show that both models are characterized by 

' the same universal scaling function. Furthermore we derive the equation of state 

^ for the pair contact process within a mean-field approach which again agrees with 

I , the mean-field equation of state of the directed percolation universality class. 

(N 
O 

' PACS numbers: 05.70.Ln, 05.50.+q, 05.65.+b 

1. Introduction 

In this paper we consider the scaling behavior of direction percolation (DP) and of 
the pair contact process (PCP) numerically in 1 + 1 dimensions and analytically 
within a mean-field approximation (see for a review [Q). Both models display a 
continuous phase transition from an active to an inactive state. In contrast to the 
unique absorbing state of DP the PCP is characterized by infinitely many absorbing 
states. Despite this difference it was observed in numerical simulations that the critical 
exponents describing the steady state scaling behavior of the PCP agree well with those 
of DP [^|, g]. This result is in agreement with a renormalization group analysis of a 
phcnomenological field theory approach of the PCP which reduces in the steady state 
to the corresponding DP equation Q] . 

In this work we consider the steady state scaling behavior of DP and the PCP in 
an external field which is conjugated to the order parameter. In particular we examine 
the scaling behavior of the order parameter (equation of state) as well as of the order 
parameter fluctuations. In contrast to previous works on the PCP we do not restrict 
our attention to the critical exponents but consider the so-called universal scaling 
functions too. In this way we are able to show that both models are characterized by 
the same universal scaling behavior, i.e., the steady state scaling behavior of the PCP 
belongs to the universality class of directed percolation. 

In the next section we briefly review the scaling behavior of DP and define 
the scaling functions. We consider then in section || for the first time the PCP in 
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a conjugated field and show that the scaling behavior is characterized by the DP 
critical exponents. The universal scaling behavior of DP and the PCP is presented 
in section [|. The obtained universal equation of state is compared to the results of a 
two loop renormalization group approach of the corresponding Langevin equation J5[ . 
Finally we derive the equation of state of the PCP within a mean-field approximation. 
Therefore we consider the PCP with particle creation at randomly selected sites. 
This random neighbor interaction suppresses long range correlations and the model is 
analytically tractable. 



2. Directed percolation 



In order to examine the scaling behavior of the 1 + 1-dimensional DP universality 
class we consider the directed site percolation process using the Domany-Kinzel 
automaton |J . It is defined on a diagonal square lattice and evolves by parallel update 
according to the following rules. A site at time t is occupied with probability p 2 (p t ) 
if both (only one) backward sites (time t — 1) are occupied. If both backward sites 
are empty a spontaneous particle creation takes place with probability p . Directed 
site percolation corresponds to the choice p x = p 2 = p and p = 0. The propagation 
probability p is the control parameter of the phase transition, i.e., below a critical 
value p c the activity ceases and the system is trapped forever in the absorbing state 
(empty lattice). On the other hand a non-zero density of (active) particles p a is found 
for p > p c . The best estimate of the critical value of directed site percolation is 
p c = 0.705489(4) [Q and we use this value throughout this work. 

The density of active sites p a is interpreted as the order parameter of the absorbing 




Figure 1. The directed percolation order parameter p a as a function of the 
particle density for zero-field (symbols) and for various values of the external field 
(h = 10~ 4 , 2 10 -4 , 5 10 -4 , 10 -3 ) (lines). The inset displays the order parameter 
fluctuations Ap a for zero field (symbols) and for various values of the external 
field h (lines). 
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Figure 2. The field dependence of the order parameter and its fluctuations 
(inset) at the critical value p c for directed percolation. The dashed lines 
correspond to the expected power-law behavior (rajjj)- 



phase transition and vanishes at the critical point according to 

Pa OC ,5/, (1) 

with Sp = (p — p c )/p c . Furthermore the order parameter fluctuations Ap a = 
l ((pI) - (Pa) 2 ) diverge as 

A Pa oc 5p-»'. (2) 

The fluctuation exponent 7' obeys the scaling relation 7' = i/ ± — 2(3, where i/ ± describes 
the divergence of the spatial correlation length at the critical point. The critical 
behavior of the order parameter is shown in figure ^ The data are obtained from 
numerical simulations of systems with periodic boundary conditions. Considering 
various system sizes L < 131072 we take care that our results are not affected by 
finite-size effects. 

Similar to equilibrium phase transitions it is possible in DP to apply an external 
field h that is conjugated to the order parameter. Being a conjugated field its has 
to destroy the absorbing phase and the corresponding linear response function has to 
diverge at the critical point i.e., 

dp* ,„s 

X = - oc. (3) 

In DP the external field is implemented (see for instance |l[]) as a spontaneous creation 
of particles (p = h > 0). Clearly the spontaneous creation of particles destroys the 
absorbing state and therefore the absorbing phase transition at all. Figure|l| shows 
how the external field results in a smoothening of the zero-field order parameter curve. 
The inset displays that the fluctuations are peaked for finite fields. Approaching the 
transition point (h — > 0) this peak becomes a divergence signalling the critical point. 
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Close to the transition point the order parameter and its fluctuations obey the 
following scaling ansatzes 

PSP, h) ~ A f. DP (*p A~^, h \~ a IP), (4) 
Ap*(5p,h) ~ A 7 ' i DP (,5pA^A CT ), (5) 

with the scaling functions f aDP and d sDP where the index sDP denotes site directed 
percolation. The scaling function of the order parameter corresponds to the equation 
of state and we recover equation (|l|) by setting Sp A -1 /' 3 = 1 at h = 0. On the other 
hand setting h \~ a IP — 1 leads to 

p a (6p = 0,h) oc h 13 ^ (6) 
at Sp = 0. Analogous we get the scaling behavior of the fluctuations at p c 

A Pa (6p = 0,h) oc ft- 7 '/- (7) 

as well as equation (||). The field dependence of the order parameter and of its 
fluctuations at p c is plotted in figure |^. 

The above scaling forms imply that curves corresponding to different values of 
the external field collapse to a single one if p a hr 13 / 17 and Ap a h 1 1° is considered as a 
function of the rescaled control parameter Sphr 1 / 17 . Using the best known estimates 
of the critical exponents (3 = 0.276486, a = 2.554216, and 7' = 0.543882 § we get 
beautiful data collapses of our numerical data (see figure ||) . 

It is worth mentioning that the validity of the scaling behavior of the equation 
of state (^) implies the required singularity of the linear response function (||) . Using 
the fact that the susceptibility is defined as the derivative of the order parameter with 
respect to the conjugated field we get 



d 

x(Sp,h) = -^p a {5p,h) 



(8) 
(9) 
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Figure 3. The scaling plot of the order parameter and its fluctuations (inset) 
for directed percolation. 
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Thus we get for the critical behavior 
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x(Sp, h) 



x(Sp, h) 



d_ 

dh 

d_ 

dh 



Pa(8p, h) 



p a (8p, h) 
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(10) 



(11) 



Sp^O 



respectively. Here the susceptibility exponent 7 is given by the scaling relation 

7 = a - (3 (12) 

which corresponds to the well known Widom equation of equilibrium phase transitions. 
Notice that in contrast to the scaling behavior of equilibrium phase transitions the 
non-equilibrium absorbing phase transition of DP is characterized by 7 ^ 7'. 

Again, the validity of the scaling ansatz of the order parameter (||) can be used 
to check that the implemented field satisfies the condition (||), i.e., to check whether 
the field is conjugated to the order parameter or not. We apply this procedure in the 
following section where we examine the scaling behavior of the PCP. 



3. Pair contact process 

The PCP as introduced by Jensen [0] is one of the simplest models with infinitely many 
absorbing states showing a continuous phase transition. At time t sites on a lattice 
of length L with periodic boundary conditions can either be occupied (ni(t) = 1) or 
empty (nj(t) = 0). Pairs of adjacent occupied sites 1, linked by an active bond, 
annihilate each other with rate p or create an offspring with rate 1 — p at either site 




Figure 4. The pair contact process order parameter p a as a function of the 
particle density for zero-field (symbols) and for various values of the external 
field from (h = 10~ 4 , 2 10" 4 , 5 10" 4 , 10" 3 ) (lines). The inset displays the order 
parameter fluctuations Ap a for zero field (symbols) and for various values of the 
external field h (lines). 
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Figure 5. The field dependence of the order parameter and its fluctuations 
(inset) at the critical value p c for the pair contact process. The dashed lines 
correspond to the expected power-law behavior (njjj)- 



i — 1 or i + 2 provided the target site is empty. The density of active bonds p a is the 
order parameter of a continuous phase transition from an active state for p < p c to 
an inactive absorbing state without particle pairs. The behavior of the PCP order 
parameter and its fluctuations are plotted in figure |]. The data are obtained from 
simulations on various system sizes L < 131072 with periodic boundary conditions. 
Our analysis reveals that the critical value is p c = 0.077093(3) which agrees with the 
value p c = 0.077090(5) fn| obtained fr om a finite-size scaling analysis of the lifetime 
distribution. 

In contrast to DP there is no unique absorbing state (empty lattice) but infinitely 
many, as any configuration with only isolated inactive particles is absorbing. Thus in 
the thermodynamic limit the system will be trapped in one of an infinite number of 
absorbing configurations for p > p c . Despite the different structure of the absorbing 
states the steady state scaling behavior of the PCP is believed to be characterized 
by the DP critical exponents f3, 7', 7 etc. On the other hand the dynamical scaling 
behavior, associated with activity spreading of a localized seed depends on the details 
of the system preparation || . 

Recently, Dickman et al. [jll] considered the PCP with an external field that 
randomly creates isolated particles. Thus the external field couples to the particle 
density but not to the order parameter itself, i.e., the external field is not conjugated 
to the order parameter. The authors observe that the external field shifts the critical 
values p c continuously and that the critical exponents are unaffected by the presence 
of the particle source. In order to investigate the PCP in a conjugated field the 
implementation of the external field of |ll| has to be modified. Several modifications 
of the external field are possible. For instance in absorbing phase transitions with 
particle conservation Jl^ | the conjugated field triggers movements of inactive particles 
which can be activated in this way E3] . 




As shown below spontaneous particle creation with rate h acts as a conjugated 
field analogous to DP. Figure |] shows how the spontaneous particle creation smoothens 
the critical zero field curves similar to the DP behavior (see figure [l]). We simulated 
the PCP at the critical value p c for various fields. The order parameter and its 
fluctuations as a function of the external field h are presented in figure pL Approaching 
the transition point, p a and Ap a scale according to the equations (g,Q) where the 
exponents {3/a and j'/a agree with the DP values. 

Furthermore we assume that the order parameter and the order parameter 
fluctuations obey analogous to DP the scaling forms 

Pa .{8p, h) ~ A fpop(«p A~ 1//3 , h A~ ff//3 ), (13) 

A Pa (5p, h) ~ A 7 ' d PCP (S P A, h A CT ), (14) 

where the distance to the critical point is now given by Sp — (p c — p)/p c - Using the 
DP values of the critical exponents j3, a and 7' we get convincing data collapses (see 
figure^). As pointed out above, the validity of the scaling ansatz ( ^3|) implies the 
singular behavior of the linear response function (g), i.e., the spontaneous particle 
creation in the PCP can be interpreted as an external field conjugated to the order 
parameter. Furthermore the data collapse confirms again that the steady state scaling 
behavior of the PCP is characterized by the DP exponents. 

4. Universal scaling behavior 

One of the most striking features of critical phenomena is the concept of universality, 
i.e., close to the critical point the scaling behavior depends only on few fundamental 
parameters whereas the interaction details of the systems do not alter the scaling 
behavior. In the case of systems with short range interactions these parameters are the 
symmetry of the order parameter and the dimensionality of space D |Q, [l5| . Classical 
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Figure 7. The universal order parameter scaling function R DP (x,l) of the 
universality class of directed percolation. The dotted line corresponds to the 
result of a two loop renormalization group analysis of the Langevin equation B 
(see text). 



examples of such universal behavior are for instance the coexistence curve of liquid- 
vapor systems Jl(| and the equation of state in ferromagnetic systems |jl7). In the 
case of absorbing phase transitions the corresponding universality hypothesis states 
that systems exhibiting a continuous phase transition to a unique absorbing state 
generally belong to the universality class of directed percolation |l9| . A different 
scaling behavior is observed only in systems with additional symmetries (such as parity 
conservation or particle conservation) or in systems with quenched randomness. 

Following the concept of universality two models belong to the same universality 
class if the critical exponents and the universal scaling functions are identical. The 
universal scaling functions R DP and Z) DP of the DP universality class can be easily 
determined by introducing non- universal metric factors Cj and di (see e.g. ^0|) for 
each scaling argument in the scaling functions (||, ||, [l3|, |l4|), i.e., 

p a (Sp, h) ~ A i? DP ( Cl Sp \- 1/p ,c 2 h \~ a IP), (15) 

A Pa (S P ,h) - A 7 ' D DP (d 1 SpX,d 2 h\ a ). (16) 

Using the normalizations i? DP (l,0) = 7? DP (0, 1) = 1 and Z) DP (1,0) = -D DP (0, 1) = 1 
the non-universal metric factors can be determined from the amplitudes of 

Pa (6p,h = 0) ~ ( Cl S P y 3 , Pa (6p = Q,h) ~ {c 2 hy/° (17) 

and 

Ap a (Sp,h = 0) ~ (d 1 5 P )-' 1 ', Ap a (Sp = 0,h) ~ {d 2 hy^l°. (18) 

Like the value of the critical point p c the non-universal metric factors may depend on 
the details of the considered system, e.g. the lattice structure, the boundary conditions, 
the used update scheme, etc. In the case of directed site percolation in the Domany- 
Kinzel automaton we have obtained the values c° DP = 2.45, & 2 DP — 0.112, d° DP = 71.6 
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Figure 8. The universal scaling function D D p(x,l) of the order parameter 
fluctuations of the universality class of directed percolation. 



and d' DP = 3984. For the PCP on a square lattice we have determined the values 
cfcp = o.665, c PCP = 0.181, d pcp = 3.21 and d PCP = 62.11. 

Analogous to the previous scaling analysis we set c 2 h \~ a IP — 1 and consider for 
both models the rescaled order parameter p a (c 2 ft.) _/3 / <T as a function of the rescaled 
control parameter c x 5p (c 2 /i)~ 1 ' <7 as well as the rescaled order parameter fluctuations 
Ap a (d 2 h) 1 / CT as a function of d 1 5p(d 2 h)~ 1 '' 7 , respectively. The corresponding data 
are presented in figure and figurejg. In both cases we get a perfect data collapse 
of the curves showing that the one-dimensional PCP steady state scaling behavior 
belongs to the universality class of directed percolation. 

Additionally to the universal scaling function R DP ^x, 1) we plot in figure^ the 
corresponding curve of an e-expansion obtained from a renormalization group analysis 
of a Langevin equation (see equation ( p5[ ) next section). Using the parametric 
representation pl| , |22| | of the absorbing phase transition Janssen et al. showed recently 
that the equation of state is given by the remarkable simple scaling function jsj 

H(9) = 6(2-9) + 0(e 3 ), (19) 

where e denotes the distance to the upper critical dimension D c — 4, i.e., e = D c — D. 
Here the scaling behavior of the quantities p a , <5p, and h is transformed to the variables 
R and 6 through the relations 

bS P = R(l-9), Pa = R?°~. (20) 

The equation of state is given by 

ah = (—) H{9) (21) 



with the metric factors a and b. The whole phase diagram is described by the 
parameter range R > and 9 e [0,2]. Using the equations (pj|,|20|,pT|) we calculated 
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Figure 9. The DP universal finite-size scaling function P and N. Both 
quantities describe the activity spreading of a localized seed (see text). 



the corresponding universal function and compare it in figure [7] with the numerically 
obtained universal scaling function R UP Or, 1). As can be seen the significant difference 
indicates that the 0(e 3 ) corrections to (fL9|) are relevant, i.e., higher orders than 0(e 2 ) 
are necessary to describe the scaling behavior of directed percolation. Of course we 
expect that this difference will decrease with increasing dimension, i.e., for e — > 0. 

Before closing this section we briefly mention that the dynamical scaling behavior 
of the PCP belongs to the DP universality class too if the spreading of a localized 
seed is considered at the so-called natural particle density ||. Examining spreading 
activity one usually considers the survival probability P a of the activity as well as 
how the number of active particles N a — Lp a increases in time. In the case of DP 
the simulations are started with a single seed on an empty lattice. For the PCP an 
absorbing state at p c is prepared to which a particle is added in order to create one 
seed (one active pair). At the critical point the following power-law behaviors are 
expected 

A a OC t\ P a (X r 5 . (22) 

Finite systems sizes limit these power-law behaviors and P a and A a obey the finite-size 
scaling ansatzes 

N a (5p = 0,L) - An(A- 1/ ^,A- 1/ez L) (23) 

P a (5p = 0,L) ~ \p(\ x / s t, X 1 ^ Sz L). (24) 

where z denotes the dynamical exponent. Analogous to the above analysis the 
universal scaling curves N and P are obtained by introducing appropriate non- 
universal metric factors. Using the values z — 1.580745, 6 = 0.313686, and 
S = 0.159464 (|] w e get convincing data collapses and the universal scaling functions 
are plotted in figure ||. 
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5. Mean field scaling behavior 

The mean-field equation of state of DP can be easily derived from the corresponding 
Langevin equation 

d t p a = 5 PPa - \pl + Kh + DV 2 p a + rj (25) 

which describes the order parameter field p a (x,t) on a mesoscopic scale (see jlj for a 
detailed discussion). Here D denotes the diffusion constant, r\ denotes a multiplicative 
noise term with the correlator 

( r]{x, t) r){x', t')) = T p a {x, t) S d (x - x!) S(t - t') (26) 

and A > 0, k > 0, and T > are certain coupling constants. Neglecting spatial 
correlations and fluctuations (D — and r\ = 0) one gets for the steady state behavior 
{dtp, = 0) 

6pp a - Xpl + Kh = (27) 
from which it is easy to derive the universal scaling function 



Pa = R DP {c^S P ,c^h) = ^ + Jdg*h + (^-) , (2S) 



where the non- universal metric factors are given by cf p = l/X and c° p = ft/A, 
respectively. For zero-field we get the solutions p a = (absorbing state) and 
p a = c° F 6p, i.e., the mean- field value of the order parameter exponent is (3 = 1. 
On the other hand 5p = leads to p a = (c° p h) 1 / 2 implying a — 2. As recently shown 
the universal scaling function (^8|) describes not only the scaling behavior of the DP 
universality class. It also occurs in the different universality class of absorbing phase 
transitions with a conserved field f23fl . 

Let us now consider the following modification of the PCP. An active bond 
produces an offspring with rate (1 — p) at an empty site selected at random. The 
rules for annihilation and action of the external field remain unchanged. This random 
neighbor interaction suppresses long range correlations and the model is therefore 
expected to be characterized by the mean-field scaling behavior. We denote the density 
of inactive bonds between an occupied and an empty site as pi. Bonds between empty 
sites are denoted as p . Normalization requires p c + p\ + p a = 1. 

Depending on the sites adjacent to the target site the number of active bonds n a , 
inactive bonds n; or empty bonds n e is changed. For instance if the adjacent sites 
are empty, for which the probability in absence of correlations is p 2 , the number of 
empty bonds decreases by two (An c = —2). On the other hand there are two new 
inactive bonds (Artj = +2). The total probability for this event is (1 — p)p a _p 2 . A list 
of all possible processes and their mean field probabilities is given in table 1. Thus 
we obtain rate equations for the expectation values iJfAn^] of the changes in active, 
inactive and empty bond numbers. These expectation values are zero in the steady 
state, i.e., 

E[An x ] = An * P{ An *) = ( 29 ) 

An, 

with x € {a, i, e}. In the case of E[An a ] we get 

E[An a ] = - 3ppl - 4pp 2 aPi - p PaP 2 + 2(1 - p)p a p 2 (30) 
+ 2(1 - p)p a pip e + 1hp\ + 2hpip e = 0, 
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Table 1. The configuration of a PCP lattice before (C) and after (C) an event. 
Only the sites left and right of those changed by particle creation (top), pair 
annihilation (middle) or particle creation due to the the external field (bottom) 
are shown. Empty sites are marked by o, and occupied sites by •. Here, An a 
denotes the change of the number of active bonds, An; the respective change of 
inactive bonds, An e that of empty bonds and P is the corresponding probability 
of the event if spatial correlations are neglected. 
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whereas we get for the for inactive bonds 

E[Am] - - 2(1 - p)p a p 2 + 2(1 - p) Pa p 2 c + 2pp\ (31) 

- 2pp a p 2 - 2hp\ + 2hp 2 = 0. 

Using this equations together with the normalization allows for calculating the order 
parameter for zero-field (h = 0) which yields the non-trivial (p a > 0) solution 



_ 8 - 3p 2 - 5 P - 2 v /2p(l - p)(3p - 4) 2 
Pa ~ V-9p + 8 ■ {6l) 

This solution is valid below the mean field critical point p c = 8/9 whereas the trivial 
solution p a = is valid for all p but unstable above p c . Expanding ( |32[ ) around the 
critical point leads to 

Pa = \tp + 0(6 P 2 ) (33) 

with Sp = (p c — p)/pc- Thus the mean field exponent of the PCP is (3 = 1 and the 
non-universal metric factor cf CP =3/8. 

In order to obtain the order parameter in presence of an external field h equations 
(^,|TJ) are solved for p\ which yields 

Ah + 4p a - Ah Pa - 4 PPa - Apl (34) 

+ {-I2p 2 pl + (2h + 2 Pa - 2hp a - 2p Pa ~ 2pl - 2ppl) 2 } 1/2 

- {(-2/i - 2p a + 2 PPa + 2p 2 a - 2pp 2 a ) 2 + ... 

... + App a {h + Pa - 2h Pa - pp a - 2p\ + hp\ + 2pp\ + pi) } 1/2 = 0. 

To obtain the field dependence of the order parameter a series expansion around h = 
at p c = 8/9 is performed which results in leading order 

3 /i, (35) 
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i.e., the mean field values of the PCP are given by c FCF = \/o/8 and a = 2. 

Finally we derive the mean-field universal scaling function R of the PCP. 
Therefore we write ( |34| ) as a function of the reduced control parameter Sp and perform 
the limits p a — > 0, — > and /i — > with the constraint that p a /vh and p a /6p are 
finite. Thus we remains in leading order with 

l8pp a - pl + ^h = 0. (36) 
Solving this equation yields 

Pa = ||«P + Qf*p) 2 = £dp« CP <^< CP ^ (37) 

Thus we have shown that both models are characterized by the same critical expo- 
nents and scaling function, i.e., similar to the 1 + 1-dimensional case the mean- field 
solution of the PCP belongs to the mean- field universality class of directed percolation. 
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